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Dust particles in orbit around a star drift toward the central star by the Poynting-Robertson effect and pile up 
by subUmation. We analytically derive the pile-up magnitude, adopting a simple model for optical cross sections. 
As a result, we find that the sublimation temperature of drifting dust particles plays the most important role in the 
pile-up rather than their optical property does. Dust particles with high sublimation temperature form a significant 
dust ring, which could be found in the vicinity of the sun through in-situ spacecraft measurements. While the 
existence of such a ring in a debris disk could not be identified in the spectral energy distribution (SED), the size 
of a dust-free zone shapes the SED. Since we analytically obtain the location and temperature of sublimation, these 
analytical formulae are useful to find such sublimation evidences. 

Key words: SubUmation - Dust - Interplanetary medium - Debris disks - Celestial mechanics. 



1. Introduction 

Refractory dust grains in orbit around a star spiral into 
the star by the Poynting-Robertson drag (hereafter P-R drag) 
and sublime in the immediate vicinity of the star. Because 
the particles lose their mass during sublimation, the ratio (3 
of radiation pressure to gravity of the star acting on each 
particle ordinarily increases. As a result, their radial-drift 
rates decrease and the particl es pile up at the outer edge 



of their subhmatio n zone (e.g.. lMukai and Yamamotoill979 



Bums et al.L 11979). Thi s is a mechanism to form a dust 



ring proposed by iBeltonI (Il966h as an accumulation of in- 
terplanetary dust grains at their sublimation zone. Ring for- 
mation of drifting dust particles is not limited to refractory 
grains around the sun but it also takes place for icy grains 
from the Edgeworth-Kuiper be lt and for dust in debris disks 
jKobavashi et all l2008l Eolol) . Therefore, dust rmg forma- 
tion due to sublimation of dust particles is a common process 
for radially drifting particles by the P-R drag. 

The orbital eccentricity and semimajor axis of a dust par- 
ticle evolve by sublimation due to an increase in its /3 ratio 
as well as by the P-R drag. We have derived the secu l ar evo - 



lution rates of the orbital elements ( Kobavashi et al. , 2009 ) 



The derived rates allow us to find an analytical solution of 
the enhancement factors for the number density and optical 
depth of dust particles due to a pile-up caused by sublima- 
tion. Our analytical solution is found to reproduce numer- 
ical simulations of the pile-up well but its applicability is 
restricted for low eccentricities of subliming dust particles. 
The analytical solution shows that the enhancement factors 
depend on dust shapes and materials as expected from previ- 
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ous numerical studies (cf. lKimura et al. . 1997 ). Although the 
solution includes physical quantities for the shapes and ma- 
terials, it does not explicitly show which quantity essentially 
determines the enhancement factors. 

The goal of this paper is to derive simplified formulae that 
explicitly indicate the dependence of dust ring formation on 
materials and structures of dust particles. In this paper, we 
adopt a simple model for the optical cross sections of fractal 
dust particles and analytically obtain not only the enhance- 
ment factors but also the location of the pile-up and subli- 
mation temperature. I n addition, we extend the model of 
Kobavashi et al. ll2009h by taking into account orbital eccen- 
tricities of subliming dust particles. 

In Section 2, we derive the sublimation temperature as a 
function of the latent heat. In Section 3, we introduce the 
characteristic radius of fractal dust and derive the sublima- 
tion distance for that dust. In Section 4, we simplify the 



formulae of enhancement factors derived bv lKobavashi et al 



(|2009) and obtain the new formulae that show explicitly the 
dependence on materials and structures of the particles. We 
provide a recipe to use our analytical formulae in Section 5, 
apply our simplified formulae to both the solar system and 
extrasolar debris disks, and discuss observational possibil- 
ities of dust sublimation in Section 6. We summarize our 
findings in Section 7. 



2. Sublimation temperature 

We consider dust particles in orbit around a central star 
with mass M*. Driven by the P-R drag, they drift inward 
until they actively sublime in the vicinity of the star. We have 
shown in lKobavashi et al. 1 20091) that the ring formation due 
to sublimation occurs only for their low orbital eccentricities 
e and obtained the secular change of semimajor axis a of the 
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particle with mass m as 



da\ 13 a dm 

'dt]^ ^^1- i3rn ~dt 



PGM, 2 
c a ' 



(1) 



where rj = — ln/3/ Inm, —dTn/dt\r=a is the mass-loss rate 
of the particle at the distance r = a, G is the gravitational 
constant, and c is the speed of light. The /3 ratio is given by 



13 



pr 



(2) 



where Cp, is the radiation pressure cross section averaged 
over the stellar radiation spectrum and is the stellar lu- 
minosity. The first and second terms on the right-hand side 
of Eq. ([T]) represent the drift rates due to sublimation and the 
P-R drag, respectively. Although we consider only the P-R 
drag from stellar radiation, the P-R drag due to stellar wind 
also transports the particles. However, the magnitude of pile- 
up, its location, and sublimation temperatu re hardly depend 



on which d rag determines their transport (iKobavashi et al. 
20Q8[l20Q9h . 



A particle generated in a dust source initially spirals to- 
ward a star by the P-R effect. As it approaches the star due 
to the P-R inward drift, its temperature rises high and it fi- 
nally starts active sublimation. The drift turns outward by 
sublimation when (3 increases with mass loss. The radial mo- 
tion of the particles becomes much slower than the P-R drift 
alone, resulting in a pile-up of the particles. Note that other 
mass-loss mechanisms such as sputtering by stellar wind and 
UV radiation are negligible during active sublimationQ 

The mass loss rate of a particle due to sublimation is given 

by 



dm 



A 



27rfcT 



Po(r) exp 



^m^H 
kT~ 



(3) 



where A is the surface area of the particle, H is the latent heat 
of sublimation, /i is the mean molecular weight of the dust 
material, m^ is the atomic mass unit, and k is the Boltzmann 
constant. Here the saturated vapor pressure at temperature T 
is expressed by Po{T) exp{— fim^H / kT) with Pq{T) being 
only weakly dependent on T. 

During active sublimation, the first term on the right-hand 
side of Eq. ([T]i increases and then (a) nearly vanishes. The 
temperature Tsub at active sublimation is approximately de- 
termined by (d) = 0. Substituting Eq. Q into Eq. ([TJ for 
(d) = 0, we have 



^sub — 



fim^H 



In 



/ 2GM, 



1-P 2TTkT, 



sub 



APo V 



(4) 

Although Eq. (|4]i is a function of a as well as Tsub, the natural 
logarithmic function on the right-hand side has little sensi- 
tivity to a and T^uh- Therefore, the sublimation temperature 



*We consider dust particles tliat can drift into their active sublimation 
zone. This is valid in the solar system, since the size decreasing timescale 
due to sputte ring is longer than the drift time due to the P-R effect 
iMukai and Sc hwehm. 1981). However, we note that icy particles may not 
come to their sublim ation zone around high ly luminous stars because of 
strong UV sputtering l lGrigorieva et al.Ll2007h . 
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Fig. 1. Subhmation temperature Tsuj, as a function of lJ.Ht; ^, where 
5 = 1 + 0.021n(Po/6.7 X 10^"' dyncm"^). The mean molecular 
weight fj., the latent heat H, and the vapor pressure in the limit of high 
temperature, Po, are listed in Table 1. The solid line indicates Eq. js). 
Circle s represent T^^h obtained from the method of .Kobayashi et alj 
120091) around the sun. 



may be approximated by 



Tsub =i 1.3 X lO^r' 



where 



H 



3.2 X IQiOergg- 



C = 1 + 0.02 In 



6.7 X lOi-^dyncm 



U] K. 

(5) 



(6) 



Here, we set m = 1.1 x lO^^^g^^ ^ 1/2, = 1/3, T^ub = 
1300 K and a = 15i?0 with the solar radius Rq in the 
argument of the logarithmic function under the assumption 
of a spherical olivine dust particle around the sun; the other 
choice of these values does not change the result significantly 
because of the slowly-varying properties of the logarithmic 
function. 

Equation (|5]) indicates that the active subhmation temper- 
ature is mainly determined by the latent heat of sublima- 
tion and mean molec ular weight of th e particles. This ex- 
plains the findings by IKobavashi et al.l ( 2008) that the tem- 
perature is insensitive to the stellar parameters, and L„. 
As a consistency chec k, we c alculate the sublimation tem- 
peratures according to IKobav ashi et al. (2009) for materials 
listed in Table 1 and compare the temperatures with Eq. (|5) 
. (see Fig.[T]i- In spite of the simplification, Eq. (|5]i is in good 
agreement with the tem perature given by the procedure of 
IKobavashi etaP ( l2009l) . 



3. Fractal dust approximation 

We introduce the characteristic radius s of a dust particle, 
which is defined as 



3 / pidV ' 



(7) 



where J dV means an integration over volume, s is the dis- 
tance from its center of mass, and pi is its interior density. 
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Table 1. Material parameters: the material density, is the mean molecular weight, H is the latent heat of sublimation, and Pq is the saturated vapor 
pressure Pv in the limit of high temperature. 





material density [g cm ^] 




H [erg g-i] 


Po [dyn cm ^] 


olivine 


3.3 


169.1 


3.21 X 10^" 


6.72 X 10" 


pyroxene 


3.3 


60.1 


9.60 X 10^" 


3.12 X 10" 


obsidian 


2.37 


67.0 


7.12 X 10^° 


1.07 X 10" 


carbon 


1.95 


12.0 


7.27 X 10" 


4.31 X lO^*' 


iron 


7.86 


55.8 


2.97 X 10^° 


5.00 X 10"^ 


pure ice 


1.0 


18 


2.83 X 10^" 


3.08 X 10^^ ^ 


dirty ice 


1.4 


18 


2.83 X 10^°'" 


2.67 X 10^^ ^ 



Note — obsidian is forme d as an igneous rock a nd may not be plausible as interplanetary dust, but is applied for a comparison with 
previous studies (e.g., Mukai and Yamamotd,ll979,) . 

and Po are obtained from the following formulae with the subli mation tempe rature of each material. 

^ Pv = 1.33 X 10* exp(-2108/T + 16.89 - 2.14 InT) dyncm"^ l lLamvlll974h 

log Pv = -2445.5646/T + 8.2312 log T - 0.016770067 + 1.20514 x 10"^T^ - 3.63227 (Pv in the cgs unit; IWashburnlll928h 



We consider that particles have a fractal structure; the mass- 
radius relation of the particles is given by to oc for a con- 
stant fractal dimension D. For the fractal dust, s — Sg\/5/3, 
where Sg is the gyration radius of the dust (IMukai et al. . 
1992h . For homogeneous spherical dust, the characteris- 



tic radius reduces to the radius of the sphere. The cross 
sections of scattering and absorption of light are approxi- 
mately described by a function of tts^ and 2tts/X, where A 
is the wavelength at the peak of light spectrum from the dust 



jMukai et aUll992h 



The smallest dust particles before active sublimation con- 
tribute most to the enha ncements of number den sity and opti- 
cal depth at the pile-up dKobavashi et alll2009l) . Small parti- 
cles produced by parent bodies in circular orbits are expelled 
by the radiation pressure if (3 > 1/2. Thus, the minimum 
characteristic radius somin of dust particles prior to active 
sublimation corresponds to (3 — 1/2. The radiation pressure 
cross section Cpr is roughly given by ttSq^-^^^ in Eq. ^ for 

Son 



> A*. Then, we have 



Son 



1.2 



Me 
AL 



P 



l.Ogcm" 



^m, (8) 



where Mq and Lq, respectively, denote the solar mass and 
luminosity. Note that p = 3TOoinin/47rSo^i„ is the effective 
density of a dust particle with the characteristic radius somin 
and mass TOomin in the following derivation. In addition, we 
discuss the application limit of our formulae in Appendix A. 
3.1 Sublimation distance 

We introduce the dimensionless parameter x. 



27rSn 



Ah ) [ l.Ogcm-3 



-1 / ^ 

/ sub 



V 1300K 



1(9) 



where Asub is the wavelength at the peak of thermal emission 
from subliming dust with temperature Tgub- We approximate 
Asub = (2898K/Tsub) /^m, which is the wavelength at the 
peak of a blackbody radiation spectrum with T^uh- 



Since we deal with dust dynamics in optically thin disks, 
the equilibrium temperature T of a dust particle at a cer- 
tain distance from a star is determined by energy balance 
among absorption of incident stellar radiation and emission 
of thermal radiation. Therefore, the relation between temper- 
ature T and distance r = a is approximately given by (e.g., 
Kobavashi et"alll2009l) 



A-KO? 



AC^<JsbT* 



(10) 



if a is much larger than the radius of the central star. Here, 
CTSB is the Stephan-Boltzmann constant and (^^(somin) and 
Cd are the absorption cross sections integrated over the stel- 
lar spectrum and the thermal emission from the dust particle, 
respectively. 

Because somin is larger than A*, the cross section 
^^(somin) is approximated by the geometrical cross section; 



C** (sOmin) 



(11) 



The cross section Cd(soniin) may be Trsgj^j^ for x ^ 1 and 
ttsqjjjJjjX for X ^ 1. We connect them in a simple form as 



C'd(sOmin) 



(12) 



When the temperature of the smallest drifting particles 
reaches Tgu b, their pile-up results in a peak on their radial 
distribution ( Kobavashi et al. [ l2009l) . With the application of 
the cross sections given by Eqs. (fTTT i, ( fT2l i to Eq. ( fTOl i, the 
sublimation distance agub at the peak is obtained as 



asub ^ 9.9 (1 



,1/2 



1/2 



1300 K 



(13) 

where Rq = 4.65 x lO^^ aU. 

Inserting x given by Eq. ^ in Eq. (fT3T l, we have Osub »: 



lV^T-^ for a; > 1 a nd a.. 



sub oc All'^h 



sub 



for X <S^ \. In 

Kobavashi et alj ( I2008I) , our simulations have shown this de- 
pendence for dirty ice under the assumption that cx AI^-^ . 
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Fig. 2. Sublimation distance a^^^^, in solar radii i?,Q as a function of 
(1 + a;-l)l/2(L*/LQ)i/2/r2^^^ where is given by Eq. (5). The 
solid line indicates Eq. jl3|. Cir cles represent a^^^^ obtained from the 
method of' Kobavashi et all (2003) around the sun. 



We coupled Eqs. (|4]i and ( fTOl l and adopted the cross sec- 
tions calculated with M ie theorjQ, and then obtained Csub 
( Kobavashi et all 1 20091) . Equatio n (IT3]) agrees we l l with 



derived from the method of Kobava shi et al. I (l2009h 



'^sub 

(see Fig. |2]i. However, Eq. (fTsT l overestimates Osub for less- 
absorbing materials (pure ice and obsidian) because our as- 
sumption of (7, = ttSqjj^jj^ is not appropriate for such materi- 
als. Nevertheless, Eq. (fT3l l is reasonably accurate for absorb- 
ing or compound dust (dirty ice). 

4. Enhancement factor 

Dust particles with mass toq in the range from rnomin to 
"lomax are mainly controlled by the P-R drag in their source 
and therefore spiral into the sublimation zone. As mentioned 
above, the smallest drifting dust with momin corresponds to 
/3 = 1/2. If the drifting timescale of dust particles due 
to the P-R drag ipR is much shorter than the timescale of 
their mutual, destructive collisions icoh the particles can get 
out of the dust source region by the P-R drag. The ratio of 
ippt to tcoi increases with mass or size. Large dust particles 
with tpR > tcoi are collisionally ground down prior to their 
inward drifts. Therefore, the largest dust momax considered 
here roughly satisfies the condition ~ tco\ at the source 
region. 

In the steady state, the number density of drifting dust 
particles is inver s ely pr oportional to the drift velocity (e.g.. 



Kobavashi et all l2009h . Since the drift velocity of dust 
particles due to the P-R drag is proportional to /?, their 
mass distribution is affected by the mass dependence -q — 
—dh\ (3 / dhim. If the differential mass distribution of the 



tWe apply the complex refractive indices of olivine from 'Huffman' f 1976) 
and Mukadand K oike 1 1990) of pyroxene from Huffman and Stann 1 1971), 
iHiioi and Takedd jT990l) . Roush et al. (1991), and Henning and Mutschkj 
<1997l) . of obsidian from Lamv (1978) and Pollack et al. (1973), of 
carb on from Manner (1987), of iron form Johnson and Christy (1974) 
and lOrdal et al. (1988) of pure and dirty ice form .Warren (.1984.) and 
iLi and Greenberd ( Il997h . 



dust source is proportional to m ^, that of dr ifting dust 
is m odulated to m^^^"^ (e.g., Moro-Martm and Malhotrii 
l2003h . Provided that successive collisions mainly produce 
dust particles in the dust source, we have h = (11 + 
3p) / (6 + 3p) for the steady state of coUisional evo lution, 
where w^/Qq cx to 



P dKobavashi and Tanaka , 2010l) . Here, 



is the specific impact energy threshold for destructive 
collisions and v is the coUisional velocity. From the hydro- 
dynamical simulations and laboratory ex periments, oc 



TO to TO 



mutations ana laboratory ex penments, ( yf, cx 
° for small dust particles dHolsapplS 1993 



Benz and Asph aug, 1999). Since p = —0.2 to for a con- 
stant V with mass, b is estimated to be 1.8-1.9. This means 
that the smallest particles contribute most to the number den- 
sity before dust particles start to actively sublime, while the 
largest particles dominate the optical depth prior to active 
sublimation. 

When the temperatures of dust particles reach Tguh, they 
start to sublime actively. Their d do not vanish perfectly, 
but they have very small |d| relative to the initial P-R drift 
velocity. The magnitude of a pile-up due to sublima tion is 
determined by the ratio of these drift rates (kobavashi et all 
20091) . Because the drift rate at the sublimation zone is 



independent of the initial mass and the initial P-R drift rate 
decreases with dust mass, the initially small dust piles up 
effectively. As a result, both the number density and the 
optical depth at the sublimation zone are determined by the 
initially smallest dust. 

The number density is a quantity that can be measured 
by in-situ spacecraft instruments, while the optical depth is a 



key fa ctor for observations by telescopes. In lKobavashi et al. 



( 20091) . we have provided enhancement factors for the num- 
ber density and the optical depth due to sublimation. Here, 
we apply the simple model for optical cross sections in 
Eqs.dTTI). (fT2] i and the properties of the fractal dust given 
by Eq. (IB.7I I. Furthermore, we take into account an in- 
crease of eccentricities from ei due to active sublimation. 
The number-density enhancement factor /at and the optical- 
depth enhancement factor fr at the sublimation zone are then 
given by (see Appendix B for the derivation) 



b — 1 — 7] 
Jn ^ —r — g{x)h{ei) + 1, 



fr 



6-1 

2 - b / TOon 
6 - 1 V TOon 



(14) 



2-6 



g{x)h{e,) + 1, (15) 



where the functions 17(3::) and h{ei) include the dependence 
on X and ei, respectively. They are given by 



2al{l + x) 
2a{l + x)+I 



(16) 



r , 1 (6-1)/?? 

h{ei) = 1 - 2 - (2/ei)^i/^ (17) 

where a — —d\n/3/d\ns = D — 2 and / = 
fj,muH / AkTsuh — 13. Since we assume that the mass 
differential number of the drifting particles is proportional 
to TOq before active sublimation, the dependence of 
fr on TOnmin/TOnmax scc n in Eq. ( fTSl l differs from that 



of iKobavashi et al. ( 2009 ^. This mass distribution is more 
realistic and consistent with that of dust particles mea - 
sured by spacecraft around the earth dOriin et al] Il985ll . 
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Fig. 3. The enhancement factors for low orbital eccentricities as a function 
of X, where the dimensionless parameter x is determined by Eqs. (9). 
Solid line represents Eq. il4l and jI5) with a use of momin = TSOmax 
and h{ei) = 1. Circles indicate the facto r numeiically calculated by 
Eq. (68) and (69) of lKobavashi et alj <2009l) around the sun for spherical 
dust listed in Table 1 . 
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10"'* 10'^ 10"^ 10"* 

Orbital Eccentricity, 

Fig. 4. Dependence of enhancement factor on ei for dirty ice, where ei 
is orbital eccentricities of dust particles at the beginning of their active 
sublimation. Solid line indicates Eq. 114) . Filled circles represent the 
results for the simulations calculated bv .Kobavashi et al. 12008.) . 



Equation ( fTTb for h{ei) is applicable for ei ranging from 
1/2^+1/ ~ 7 X 10"*^ to 1/2/ ~ 0.05. Dust particles do 
not pile up for ei > 1/2/ and hence we give h{ei) — 
for ei > 1/2/ dKobavashi et all [20091) . In addition, 
h{ei) = 1 for ei < 1/2^"*"^/, while drifting dust particles 
hardly reach such small eccentricities (ei < 1/2*^+^/ ^ 
10~^) because their eccentricities naturally rise as high as 
the ratio of the Keplarian velocity to the speed of light 
[- 10"^(a/l AU)-i/2(^//^/^/q)1/2] by jjjg p.R gffggj^ 

In Fig. [3] we compare the simplified formulae given by 
Eqs. (fT4l) . (fTSl ) with the enhancement factors rigorously cal- 
culated by the formulae of Kobayashi et al. (2QQ2). The x 
dependence of the enhancement factors given by Eqs. (O, 
( fTSl ) is shown in the function g{x), which is an increasing 
function ranging from 2a (x = 0) to I (x — oo). Equa- 
tions (O, ( fTSl ) briefly explain the tendency of the enhance- 
ment factors; the materials with high x produce high en- 
hancement factors. 

Because x is proportional to Tguh/p (see Eq. (|9]l), the 
enhancement factors increase with Tguh/p- Thus, materi- 
als with high sublimation temperature tend to pile up suf- 
ficiently. In addition, fluffy dust particles wit h D 2 ca n- 



not effectively pile up even though p is low (iKimura et al 
[r997). This is explained by the low a = D — 2m the func- 
tion g{x). Particles with D = 3 produce the highest en- 
hancement factors. In spite of D = 3, compact particles are 
not the best for the pile-up due to a high density. Dust parti- 
cles composed by ballistic particle-cluster aggregation have 
/) ~ 3 but low effective densities relative to compact ones. 
Therefore, such porous particles with /) ~ 3 may produce 
high enhancement factors due to large x resulting from their 
low densities. In addition, high x around a luminous star 
brings the enhancement factors to increase with stellar lumi- 
nosity, which is shown for di rty ice, obsidian, and carbon in 



Kobayashi et al. V2008','2009^. 

In [Koba yashi et al. (2008) , we show the eccentricity de- 
pendence of enhancement factors from our simulations. The 



dependence is explained by h{ei) in the simplified formu- 
lae (see Fig. |4|i. Dust particles can pile up sufficiently for 
ei < 10~'^ because of h{ei) ~ 1. Otherwise, the enhance- 
ment factors decrease with ei. For ei > 0.05, the sublima- 
tion ring is not expected. 

5. Recipe 

We briefly show a recipe to obtain the sublimation tem- 
perature Tsub, its distance Csub, and the enhancement factors 
Jnt It- At first, the sublimation temperature Tgub is avail- 
able from Eq. (|5]l adopting the material properties /i, H, 
and Po listed in Table 1. Then, we calculate the dimen- 
sionless parameter x through Eq. (|9]i, applying the stellar 
luminosity and mass of interest and the bulk density listed 
in Table 1 for compact spherical dust. Note that we should 
adopt a lower density for porous particles, taking into ac- 
count their porosity. Inserting x in Eq. (fT3b . we derive the 
sublimation distance Ogub- We further need orbital eccen- 
tricities ei of dust particles at the beginning of active sub- 
limation to calculate the enhancement factors, Jn and fr- 
The dust particles resulting from collisions have eccentric- 
ities eo ~ /? at the distance ao of the dust production re- 
gion. Since particles with the highest (3 contribute most to 
a sublimation ring, we estimate eo ^ 0.5. Because eccen- 
tricities are dumped by the P-R drag, we can calculate ei 



from the relation asubS 



-4/5, 



1 



-4/5/1 



dWvatt and WhippleLll95d) . Inserting x and ei to Eqs. (IT4l) 
and ( fTsT l. we obtain and fr- 

6. Discussion 

The asteroid belt and the Edgeworth-Kuiper belt (EKB) 
are possible dust sources in the solar system. A dust counter 
on board spacecraft can measure the number density of dust 
particles. The sublimation of icy dust occurs at Osub = 
20 AU given by Eq. ( fT3] ). Icy particles reaching the sub- 
limation zone from the EKB still have high eccentricities 
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ei 



> 



0.1 ("Kobavashi et al., 2008). Therefore, a substan 



tial sublimation ring is unexpected because of — 1 for 
h{ei) = 0. Since the number density of dust particles de- 
creases inside the sublimation zone, only a bump in the ra- 
dial profile of the num ber density appears around Csub (see 
Kobavashiet"alll2010l) . In contrast, dust particles composed 
of rocky, refractory materials actively sublime at several so- 
lar radii from the sun. Therefore, orbital eccentricities of 
dust particles coming from the asteroid belt drop to ^ 0.01 
around the sublimation distance due to the P-R drag. Since 
we have /i(O.Ol) ~ 0.3 from Eq. ^ for 6 = 11/6 and 
D = 3, fN ^ 1.3-3.0 is obtained form Eq. (0 for a; > 1. 
Therefore, future in-situ measurements of dust could find 
such a sublimation ring of refractory dust particles originat- 
ing from the asteroid belt, but not a ring of icy dust particles 
from the EKB. 

A dust ring was observed around 47?© from the sun in the 
period 1966-1983, althoug h it was not detected in the 1990s 
( Kimura and Manil 1998, for a review). The optical depth 
is measured by dust emission observations. The enhance- 
ment factor fr of the observed ring is e stimated to be 2-3 



dMacOueeni Il968l: iMizutani et al.L Il984i) . The mass distri- 
bution of drifting dust particles for b = 11/6 is consistent 
with the measurement of dust particles with masses ranging 
from TOomi 



10- 



^ g by spacecraft 
Using that, we 



10' 

around the earth orbit (iGriin et al.L Il985l) 
estimate < 1.1 from Eq. ( fTST i for ei ~ 0.01. Thus, the 
enhancement by sublimation cannot account for the observed 
dust ring. However, this low fr is mainly caused by the mass 
range of drifting dust particles. If momax decreased during 
the transport of dust particles from the earth's orbit to the 
sublimation zone, higher fr could be expected. For example, 
the largest particles in the mass distribution become smaller 
by the sputtering from the solar wind. The sputtering may 
decrease ei as well as momax- Small particles with high ec- 
centricities from the dust source are ground down by sputter- 
ing and blown out by the radiation pressure before reaching 
the sublimation zone, while large particles with low eccen- 
tricities gradually become small by sputtering without the 
increase of their eccentricities and drift into the sublimation 
zone. If the ratio of tpR to the timescale of decreasing size 
due to sputtering ranges in 0.1-0.7, sublimation could form 
such a bright ring becau se of small mpma x ^Pd e i. Indeed, 
the ratio derived by Mu kai and Schwehm ( 198ll) is consis- 
tent with the condition for the formation of a sublimation 
ring. That may be a clue to explain the observed ring. 

Debris disks found around main sequence stars would 
be formed through collisional fragmentation in narrow 
planetesimal belts, which may resemble the asteroids and 
Edgeworth-Kuiper belts in the solar system. In young de- 
bris disks, fragments produced by successive collisions are 
removed from the disk by radiation pressure. We call such 
a disk a collision-dominated disk. Once the amount of bod- 
ies has significantly been decreased through this process, the 
P-R drag becomes the main removal process of fragments. 
Such a disk is referred to as a drag-dominated disk. We have 
investigated the dust ring formation in drag-dominated disks. 
To observe a sublimation ring requires a high enhancement 
factor for the optical depth. As shown in Eq. ( fTST i, a 
small ratio of momax to momin yields high fr. The condi- 



tion of TOomax Womin IS cxpcctcd to form a bright ring. 
Since the drift time due to the P-R drag is comparable to 
the collisional time for bodies with momax, the condition of 
"^Omax ^ ?TiOmin IS achieved in transition from a collision- 
dominated disk to a drag-dominated one. 

A significant sublimation ring consisting of icy particles 
is not expected in a debris disk due to high eccentricities if a 
planetesimal belt as a dust source is located within a few hun- 
dreds AU, similar to the solar system. On the contrary, dust 
particles composed of refractory materials have ei ~ 0.01 
or smaller if a planetesimal belt is around the distance of 
the asteroid belt or further outside. Furthermore, refractory 
dust particles have high sublimation temperatures and hence 
produce a higher enhancement factor Recently, inner de- 
bris disks of refractory grains have been observed through 
interferomet ry around Vega, r Get, C Aql, and Leo, and 



Akeson et al 



Formalhaut ( Absil et all 120061 120081: iDi Folco et al.L 12007 



2009h . Such inner debris disks may have no 



table sublimation rings. 

To check the observability of a sublimation ring in the 
spectral energy distribution (SED) of thermal emission ex- 
pected from a disk around Vega located at the distance of 
7.6 pc from the earth, we take our formulae with Af* = 
2.1 Mq and L^. = 59 Lq for compact spherical olivine par- 
ticles. We obtain T^uh = 1300 K, Usuh = 0.35 AU and 
fr = 3.1, where we adopt momax = "lomin and h{ei) = 1 
in Eq. ( fT5b . The smallest radius somin = 10 /xm is much 
larger than the peak wavelength of thermal emission ( Agub — 
2.3 /im) and hence we simply treat dust particles as black- 
bodies to calculate the SED. The optical depth of dust 
particles drifting by the P-R drag without sublimation from 
the outer edge Oout to the inner one is given by a con- 
stant To, where equals Csub- The total optical depth t of 
the disk can be set as 



^ f To for flin < r < Oout, 

I otherwise, 
Toifr - l)(asub + Sasnh " r)/Sasuh 

for flsub <r < Osub + Sasnh, 
otherwise, 

(18) 

where To denotes the optical depth of dust particles form- 
ing a dust ring with width (5aaub- The width of the ring is 
roughly given by Susuh = 0.05 Osub ("Kobav ashi et al. , 20081 
2009). We set Oout — l-5asub- Note that the other choice 
of Oout does not change our result drastically. Choosing 
To = 8 X 10 ~'^, we c an reproduce the flux density ~ 8.7 Jy 



measured bv I Absil et al. (2006) with the interferometry at a 
wavelength of 2.1 /zmH Figure |5]depicts that the sublimation 



-fNote that tq depends on the collision and drift timescales, tc, <pr, in 
a planetesimal disk. Drag-dominated disks should satisfy the condition 
tc ^ tpR', TO < (1 + 7)i'k/c with the Keplarian velocity 1)]^ and the 
ratio 7 of the P-R drag force due to the stellar wind to that due to the 
stellar radiation. The value of tq applied for fitting is much larger than 
that for the drag-dominated disk (to < 3 x 10~*), if we only consider 
the P-R drag by the stellar radiation. However, the mass loss rate of Vega 
is estimated to be l ess than 3.4 X 10"^'' Af0yr~^ from radio-continuum 
observations (HoU is et allll985 5). For the upper limit of the mass loss rate, 
drag-dominated disks can have the optical depth to < 7 X 10~^ because of 
the P-R drag due to the stellar wind (7 ~ 300). Thus the disk around Vega 
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Fig. 5. The spectral density distribution of the dislc around Vega. The 
solid line indicates the flux density from a disk with the sublimation ring, 
choosing tq = 8 X 10"'^ to fit the observational data at 2.1 /.tm. Dot- 
ted lines represent that without a sublimation ring (/t = 1 resulting in 
To = 0) for TO = 1 X 10~^. To check the a^^-ah dependence, we set 
ciiii = O.Gagui, for tq = 4 X 10"'' in the case without the sublimation 
ring (da shed lines). Circle corresponds to the interferometric measure- 
ment by Absil et al. 1 2006) and diamonds indicate the photometric data 
jAbsil et al..2006- and reference therein). 



ring does not bring about a noticeable spectral feature in the 
SED, while the SED strongly depends on a^uh value. The 
flux density from the disk diminishes with decreasing wave- 
length for the wavelength smaller than Asub because of the 
absence of dust particles with temperature higher than T^nh 
due to sublimation. The result for a-m — 0.6asub — 0.21 AU 
is shown to better agree with the observational data. Since 
flsub ^ 0.22 AU (Tsub =i 1700 K) for pyroxene, the Vega 
disk seems to be abundant in pyroxene compared to olivine. 
That could be recognized as an evidence of sublimation un- 
less the light scattering of dust particles exceed their thermal 
emissionl^ 

7. Summary 

1. We provide formulae of enhancement factors for the 
number density and the optical depth due to a pile-up 
caused by sublimation and its location and sublimation 
temperature, applying a simple model for optical cross 
sections of fractal dust particles. 

2. High sublimation temperatures result in substantial en- 
hancement factors, though the pile-up is insensitive to 
the optical properties of dust particles. 

3. If we adopt the mass distribution measured around the 
earth, the enhancement factor for the optical depth near 
the sun is smaller than 1.1. Therefore, the enhance- 
ment cannot explain the solar dust ring detected in the 
epoch of 1966-1983, unless the largest particles were 
destroyed by sputtering. 



may be a drag-dominated disk. 



§The scattering of li ght from the disk a round Vega is negligible around the 
wavelength ~ 1 um lAbsil et al.Ll200d) . 



4. The number-density enhancement factor is expected to 
be 1 .4-3 in the vicinity of the sun. Therefore, a subli- 
mation ring could be found by in-situ measurements by 
spacecraft around several solar radii from the sun. 

5. Sublimation removes dust particles with temperatures 
higher than their sublimation temperature. In the spec- 
tral energy distribution, the flux density from a disk re- 
duces with decreasing wavelength, if the wavelength is 
shorter than the peak one of the blackbody spectrum 
with the sublimation temperature. That could be seen 
as a sublimation evidence. However, it is difficult to 
find signs from a sublimation ring in the spectral energy 
distribution. 

Appendix A. Application limit 

The value of (3 increases with decreasing radius as long 
as the radius fulfills the condition s > A*. For s < A*, 
however, it decreases with decreasing radiusjj] Hence, /3 has 
a maximum value at s ^ A* = (2898 K/T*) fim. From 
Eq. (|2]i, the maximum value of j3 is approximately given by 



max -0.98 ( ^ 



Mq 



5 X lO^Ky/ Vl-Ogcm 



P 



(A.l) 



where the radiation pressure efficiency averaged over the 



stellar radiation spectrum Cpr 
define the effective density by 



P - 



3m 



for s > A,. Here we 



(A.2) 



with the use of the dust mass m and the characteristic radius 
s. Note that p depends on s in general; p is constant for 
(e.g., a compact sphere), whereas p is proportional 



to for a fracta l aggregate with fractal dimension D (see 



Mukai et al.Lll992L for the relation) 



When dust particles are produced by successive collisions 
between large bodies, their largest /3 does not exceed 1/2 
because the dust with /? > 1/2 cannot resist against strong 
radiation pressure. These particles then drift into the subli- 
mation zone. The smallest drifting dust particles that have 
the largest /3 contribute most to the pile-up caused by sub- 
limation. Because the enhancement factors of the number 
density and the optical depth are proportional to the largest 
/? value that drifting dust particles attain, we may expect in- 
sufficient pile-ups of the particles if /Smax < 1/2- We thus 
derive simplified formulae for characteristics of a dust ring 
from the assumption of /Smax > 1/2. Because the variation 
of is much larger than that of for main sequence stars, 
this assumption is translated into the condition for stellar lu- 
minosity given by 



> 0.5 



5 X 103K 



P 



l.Og cm 



if Cpr — TTs^ for s > A*. 



(A.3) 



I ^ Not e that /3 is independent of s for much smaller particles l lGustafsonl 
1 1994). 
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Appendix B. Derivation of enhancement factors 

We assume that the mass differential number of drifting 
dust particles is proportional to ttiq for the drifting parti- 
cles with mass mo. In addition, we adopt /3(m) oc m^^', and 



S{m) 



where S is the geometrical cross sectio n of 



a dust particle with mass m. In lKobavashi et alj | |2009[) . we 
derived the number-density enhancement factor /at and the 
optical-depth enhancement factor at the peak as 



In 

fr 



??("l0min)5m(?lub,TO0min) ^-2(2/!: 



2/2 j 

y2) 



where 



^1(2/1,2/2) 



^■2(2/1,2/2) 



/r y-'dy 



l,(B.l) 
1,(B.2) 



(B.3) 



(B.4) 



As we 



2/1 = m-init,max/"10min, and y2 TOomax/"^Omln■ 

will describe below, dust particles with niQ < mi„it,max 
drifting into the sublimation zone can contribute to the en- 
hancement factors. 

Here g^n is a function of the optical properties of dust 
particles with m = momin at Tsub, namely, given by a 
function of x. The function gm is defined as 



i(r,m) 



4+ct 



^ d In Py 
k dlnT 



n 

1-p 



where 



CT 



dlnm ' 



Cd 



dlnr) 
d In m 



91nCd 



' d In Py 
V dlnT 



(l2009l) derive the relation between e and m during the active 
sublimation as 



l-/?(mi) 
1 - /3(to) 



ei, 



where 



1 



4 + CT V d\nT 



(B.9) 



(B.IO) 



with mi is the mass starting the active sublimation^ We 
approximate k = I because / ^ 1. Equation ( IB. 91 ) in- 
dicates that e substantially changes for (3^1. Since the 
mass loss is negligible outside the sublimation zone, we ap- 
proximate nil c:^ fno- Partic les can pile u p as long as e < 
2kTsub/tJ-muH ^ 1/21 (Kob avashi et al.L l2009). Substitut- 
ing e = 1/2/, /3 = 1/2, and = ;S(minit,max/"l0min)"''' 

into Eq. (IB.9I I. we have 



^init,max 



(2/ei) 



2(2/ei)i// _ 1 



momin- (B.U) 

< TOomax and Ci > 



Equation dB.l II ) is valid for mjnit.max 
l//2^+^. We should set mi„it,max — "T-Omax instead of 
Eq. ( IB. lib for mi„it,max > "lomax or ei < l//2^+^ Be- 
cause y2 > 1, and b = 1.8-1.9 in Eqs (IB. 3b and ( IB.4b . hi 
and /i2 are given by 



2 J 4+CT 

(B.5) 



91ni 



dlnA 
d In m 



(B.6) 



According to our simple model in Sections|2]and|3] a particle 
with Tsub and momin has 

= 2/D, Cd = 2/D + l/{x + 1)D, Ca = Cs = 2/D, 
Tj = -d\nP/dlnm= {D -2)/D, /3 = 1/2. (B.7) 



hiiyi; 
^2(2/1- 



2/2) = 



2/2 j 



b-1 
2 — b / mo 



/i(ei), 



1 V"^On 



2-6 



/i(ei 



(B.12) 



(B.13) 



where h is defined as Eq. ( fTTb . 

Substituting Eqs. (lOl i. ( IRT2] i. and dBJJb into Eqs. JbTT) 
and (IB. 2b . we have the enhancement factors in Eqs. ( fT4l i and 
dBJ. 
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gm{Tsuh: "T^Omin) 



2ID{l + x) 
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